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PREFACE 


The  overall  objective  of  the  project  was  to  measure  and  analyze  the 
global  ice-hull  interaction  forces  of  a  ship  ramming  in  multi-year  ice.  Spe¬ 
cific  objectives  were  to  obtain  a  larger  data  base  of  ice  loads  for  the 
development  of  analytical  models  describing  ship-ice  interaction  and  analyze 
the  effect  of  ship  displacement  and  bow  shape  on  global  ice  loads  using  the 
Canadian  experience  gained  from  the  1983  KIGORIAK  and  ROBERT  LEMEUR  Impact 
Tests  [13]. 

An  instrumentation  system  was  developed  for  measurement  of  global  ice 
impact  loads  onboard  the  USCGC  POLAR  SEA.  From  September  29,  1985,  to 
October  12,  1985,  the  POLAR  SEA  conducted  ice  impact  tests  on  heavily  ridged 
ice  features  in  the  Alaskan  portion  of  the  Beaufort  Sea.  8ending  strain  gage 
measurements  were  used  to  estimate  the  longitudinal  bending  moment  distribu¬ 
tion  of  the  POLAR  SEA  during  an  impact  with  an  ice  pressure  ridge.  Compres¬ 
sive  strains  along  the  stem  and  ship  acceleration  and  velocity  measurements 
were  also  recorded.  Hull  stress  and  impact  force  time-histories  were  calcu¬ 
lated  along  with  the  longitudinal  bending  and  shear  distributions  during  ice 
impact  events.  The  results  indicated  that  the  methodology  used  in  estimating 
the  impact  force  provided  a  greater  understanding  of  the  ship-ice  interaction 
process . 


One  of  the  main  conclusions  from  this  study  was  that  the  ice  force  on 
the  bow  of  the  POLAR  SEA  should  not  be  treated  as  a  localized  load.  The  load 
is  spread  over  a  large  area  of  the  bow  from  the  stem  hook  forward.  The 
longitudinal  location  of  the  maximum  bending  stress  was  just  forward  of  the 
super-structure.  The  dominant  response  of  the  vessel  was  first  mode  vibration 
at  about  3  Hz.  The  loading  rate  was  measured  to  be  as  high  as  50Q0  LT/sec  (50 
MN/sec),  wnicn  is  much  lower  than  the  15,000  LT/sec  (150  MN/sec)  loading  rate 
noted  on  the  KIGORIAK  [13].  The  maximum  vertical  bow  force  observed  during 
ramming  was  2506  LT  (24.97  MN)  with  a  maximum  measured  bending  stress  of  6,078 
psi  (41.91  MPa).  This  is  well  below  the  45,600  psi  (310  MPa)  yield  strength 
for  the  hul 1  materi al . 


It  is  recommended  that  further  analysis  be  done  to  determine  the  shape 
of  the  pressure  distribution  during  each  impact  and  to  estimate  the  transverse 
force  from  unsymmetric  rams.  Both  can  be  determined  by  additional  analysis  of 
the  measured  data.  It  is  further  recommended  that  additional  multiyear  ice 
data  be  collected  with  the  POLAR  Class.  Additional  data  is  required  to  better 
understand  the  relationship  between  vertical  bow  force  and  impact  velocity. 
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1.  INTRODUCTION 


The  USCG  POLAR  Class  winter  deployments  sponsored  by  the  Maritime 
Administration  (MarAd)  have  provided  a  platform  from  which  to  gather  environ¬ 
mental,  trafficabi 1 ity ,  and  ship  performance  data.  For  this  phase  of  the  pro¬ 
gram,  a  deployment  of  the  POLAR  SEA  in  September  and  October  of  1985,  the  Ship 
Structure  Committee,  the  American  8ureau  of  Shipping  and  the  Maritime  Adminis¬ 
tration  in  conjunction  with  the  Canadian  Ministry  of  Transport  sponsored  a 
program  to  collect  global  ice  loads.  Additional  support  for  this  program  as 
well  as  the  concurrent  environment  data  collection  program  was  provided  by 
Newport  News  Shipbuilding,  Bethlehem  Steel  and  five  members  of  the  Alaskan  Oil 
and  Gas  Association  (AOGA),  Arco,  Chevron,  Exxon,  Mobil  and  Shell. 

The  ultimate  objective  of  this  jointly  funded  research  is  to  develop 
ice  load  criteria  for  the  future  design  of  ships.  Specifically,  the  objective 
of  this  study  was  to  measure  the  total  load  that  ice  exerts  on  the  hull  of  the 
vessel  when  it  rams  large  ice  features.  Other  objectives  included  increasing 
the  data  base  of  ice  loads  for  the  development  of  analytical  models  describing 
the  ship-ice  interaction  and  for  understanding  the  effect  of  ship  displacement 
and  bow  shape  upon  the  global  ice  loads  by  comparison  with  other  available 
data. 


The  "global  ice  load"  is  defined  as  the  net  resultant  of  the  ice  loads 
generated  at  the  many  local  contact  areas  around  the  bow  during  impact.  These 
loads  may  generate  significant  bending  moments  in  the  hull  girder,  which  may 
affect  the  structural  integrity  of  icebreaking  ships.  This  in  turn  has 
implications  on  the  design  of  icebreaking  vessels  and  the  type  of  design 
criteria  to  be  developed  [1], 

Since  the  start  of  commercial  oil  development  in  the  Arctic  a  number  of 
analytical  models  describing  ship-ice  interaction  have  been  developed  using  a 
rigid  body  idealization,  flexible  beam  elements,  and  three  dimensional  finite 
element  models  [2, 3, 4, 5].  Full-scale  impact  tests  have  also  been  con¬ 
ducted  on  the  icebreaking  vessels  M.V.  CANMAR  KIGORIAK  [6,7],  M.V.  ROBERT 
LEMEUR  [7,9],  M.V.  ARCTIC  [10]  and  now  the  USCGC  POLAR  SEA.  General 
discussions  of  these  tests  can  be  found  in  references  8  and  11.  Physical 
modelling  of  the  ship-ice  impact  interaction  for  the  M.V.  ARCTIC  has  also  been 
carried  out  by  ARCTEC  CANADA  for  the  Canadian  Coast  Guard  [12]  and  the  Tech¬ 
nical  Research  Center  of  Finland  under  a  joint  research  program.  The  focus  of 
all  the  work  has  been  to  provide  a  sound  technical  basis  for  further  develop¬ 
ment  of  ice  load  design  criteria  to  accommodate  the  technical  and  regulatory 
requirements  of  expanding  maritime  operations  in  the  Arctic. 

The  work  presented  here  was  carried  out  onboard  the  USCGC  POLAR  SEA  in 
the  Alaskan  Beaufort  Sea  between  September  19  and  October  13,  1985.  This 
report  describes  the  way  the  global  ice  loads  were  collected  as  well  as  a 
presentation  and  analysis  of  the  collected  data.  Figure  1  shows  the  principal 
areas  of  operation  during  the  deployment. 


2 


153°  157°  156°  155°  154°  153° 


i  i  i 


Figure  1 

DAILY  LOCATION  OF  POLAR  SEA 
AT OOOO  LOCAL  TIME 
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2. 


DESCRIPTION  OF  THE  OATA  ACQUISITION  SYSTEM 


An  outline  or  the  instrumentation  system  developed  for  this  project  is 
illustrated  in  Figure  2.  The  present  system  was  adapted  from  the  system 
originally  developed  by  Canadian  Marine  Drilling  Ltd.  (Canmar)  for  the  1983 
M.V.  KIGORIAK  and  M.V.  ROBERT  LEMEUR  full-scale  impact  tests  [7J.  There  are 
several  fundamental  differences  between  these  two  systems.  The  approach 
used  for  Polar  Sea  measured  the  longitudinal  bending  strain  distribution  span¬ 
ning  the  location  of  the  ice  load,  whereas  the  Canmar  system  measured  the 
shear  strain  in  sections  near  the  location  of  the  ice  load  and  the  bending 
strain  distribution  aft  of  the  location  of  the  ice  force.  The  POLAR  SEA 
method  required  fewer  strain  gages  at  each  location  and  therefore  allowed 
more  frames  along  the  ship  to  be  instrumented.  The  result  was  an  excellent 
definition  of  the  longitudinal  bending  and  shear  distribution  spanning  the 
location  of  the  load  because  of  the  larger  number  of  instrumented  frames. 
Secondly,  in  these  tests  the  actual  longitudinal  location  of  the  load  was 
measured  from  compressive  strain  gages  along  the  centerline  bulkhead,  while 
the  Canmar  system  had  to  infer  the  location  from  other  data.  Additional 

details  of  the  POLAR  SEA  system  as  well  as  an  itemized  channel  description 

are  presented  in  Appendix  A. 

To  estimate  the  vertical  ice  force  on  the  bow  during  an  impact  with  a 
heavy  ice  feature,  the  shear  force  around  the  location  of  the  load  must  be 
well  defined.  Figure  3  gives  some  idealized  Shear  and  bending  moment  diagrams 
for  an  icebreaker  ramming  into  an  ice  feature.  As  the  lower  figure  indicates, 
the  shear  force  changes  from  negative  to  positive  over  a  relatively  short 
distance  around  the  location  of  the  load.  Since  the  shear  force  is  the 
negative  of  the  slope  of  the  bending  moment  diagram,  then  the  bending  moment 
must  be  well  defined  over  this  same  region  in  order  to  obtain  an  accurate 
estimate  of  the  shear  force.  With  this  in  mind,  the  majority  of  the  frames 
instrumented  for  bending  were  concentrated  near  the  anticipated  location  of 
the  ice  force.  Figure  4  shows  the  location  of  these  gages. 

The  bending  gages  along  the  01  Deck  and  on  the  3rd  Deck  at  frame  55 

were  placed  parallel  to  the  side  shell  in  pairs  along  opposite  sides  of  the 
ship.  Measurements  taken  from  these  gages  were  later  transformed  into  the 
strain  parallel  to  the  centerline.  (See  Appendix  B  for  the  details  involved 
in  any  of  these  conversions  and  computations.)  In  the  calculation  of  the 
longitudinal  strain  due  to  the  ice  force,  the  data  from  each  pair  were 
averaged  together  to  exclude  any  torsional  strain.  Another  advantage  of  this 
gage  pair  arrangement  was  the  ability  to  observe  the  symmetry,  or  the  lack 
thereof,  in  the  ice  loading  during  a  ram. 

The  bending  gages  were  placed  on  at  least  two  levels  for  every  location 
forward  of  frame  86.  This  arrangement  allowed  the  computation  of  the  bending 
moment  based  upon  a  stress  couple  that  was  a  known  distance  apart.  It  had 
the  further  advantage  of  eliminating  the  longitudinal  stress,  and  therefore 
force,  from  the  bending  moment  calculation. 
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Figure  3 

ESTIMATED  SHEAR  AND  BENDING  MOMENT  DIAGRAMS 
FOR  THE  USCGC  POLAR  SEA 
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CENTER  LINE  PROFILE 


LOCATION  OF  LONGITUDINAL  COMPRESSION  GAGES 
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Compression  gages  for  estimating  the  location  of  the  load  were 
installed  from  cant  frame  14  to  cant  frame  38.  They  were  placed  on  the  cen¬ 
terline  bulkhead,  just  forward  of  the  cant  frames  and  12  to  18  inches  above 
and  perpendicular  to  the  stem  bar.  Figure  5  shows  these  locations.  The  spac¬ 
ing  of  these  gages  enabled  an  accurate  estimate  of  the  center  of  the  ice 
force  to  be  made  since  the  load  could  be  "sensed"  every  32  in  (80cm)  along  the 
stem.  The  gage  distance  away  from  the  stem  bar  was  selected  from  work  done  on 
the  placement  of  gages  for  measurement  of  local  loads  [15]  to  avoid  the  possi¬ 
bility  of  "dead  spots"  between  the  gages.  This  system  also  provided  an  esti¬ 
mate  of  the  impact  speed  as  the  peak  ice  force  moved  along  the  stem  bar. 

Transverse  compression  gages  were  installed  at  cant  frame  35  as  shown 
in  Figure  6.  The  output  from  these  gages  was  intended  for  use  in  estimating 
the  width  of  the  ice  contact  area  and  the  distribution  of  pressure  across  the 
bow  of  the  POLAR  SEA.  These  gages  were  added  knowing  that  the  analysis  of  the 
data  could  not  be  completed  within  the  intitial  scope  of  work.  It  is  hoped 
that  this  data  can  be  reduced  in  the  future. 

In  addition  to  the  strain  gages,  three  uniaxial  accelerometers  arranged 
in  a  triaxial  array  and  oriented  to  the  ship's  principal  axes  were  located  in 
the  bow  area  as  shown  in  Figure  4.  The  output  from  the  yaw  accelerometer  was 
used  to  determine  if  a  ram  was  symmetric.  The  accelerometer  readings  could 
also  be  used  in  future  analysis  to  provide  an  estimate  of  the  inertial  forces 
forward  of  the  ice  force  and  help  assess  the  relative  importance  of  the 
longitudinal ,  transverse  and  vertical  ice  forces  acting  on  the  vessel.  The 
POLAR  SEA  was  also  equipped  with  a  doppler  microwave  speed  log.  This  radar 
was  mounted  in  the  waist  of  the  vessel  and  oriented  forward  to  provide  an 
estimate  of  the  impact  velocity.  The  specifications  and  locations  of  all  the 
transducers  used  for  the  onboard  instrumentation  are  described  in  Appendix  A. 

The  required  sampling  frequency  for  measurement  of  the  strain  response 
of  the  POLAR  SEA  was  selected  based  on  the  rate  of  loading  and  the  vibrational 
frequency  of  the  ship.  Previous  experiments  have  indicated  that  the  dominant 
vibrational  frequency  was  approximately  3  Hz  [16].  The  predicted  rise  time  of 
the  ice  force  was  used  to  estimate  the  rate  of  loading.  In  this  case, 
previous  full-scale  measurements  indicated  rise  times  to  be  as  fast  as  0.1 
seconds  [17].  If  a  quarter  sine  wave  is  assumed  for  the  rise  in  strain,  a 
corresponding  maximum  frequency  of  interest  of  2.5  Hz  results  (period  of  0.4 
seconds).  A  low-pass  filter  frequency  of  10  Hz  was  selected  such  that  it  was 
well  above  all  the  frequencies  of  interest.  The  minimum  digital  sampling 
frequency  would  then  be  32  Hz  to  ensure  a  unique  10  Hz  sine  wave.  This  is 
exactly  the  system  that  was  used  in  the  local  loads  measurement  program  [15]. 
In  this  case  a  more  sophisticated  data  acquisition  system  allowed  an  increase 
in  sampling  frequency  over  the  local  loads  system,  so  100  Hz  was  selected  to 
provide  at  least  10  samples  during  the  strain  rise  time.  Data  was  sampled  for 
25  seconds  which  was  determined  by  the  size  of  the  storage  medium.  To 
increase  sampling  frequency  beyond  100  Hz,  the  length  of  recording  would  have 
to  be  shorter  or  a  larger  storage  medium  could  be  used.  The  chosen  values 
were  a  reasonable  compromise  for  the  intended  measurements. 
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Figure  6 

LOCATION  OF  TRANSVERSE  COMPRESSION  GAGES 
ON  CANT  FRAME  35 
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3. 


DESCRIPTION  OF  THE  TEST  PROGRAM 


Global  load  data  were  gathered  on  the  summer  deployment  of  POLAR  SEA  to 
the  Beaufort  Sea.  The  program  was  completed  in  conjunction  with  an 
environmental  data  collection  effort  involving  extensive  on-ice  surveying  of 
surface  and  sub-surface  topography.  Measurement  techniques  included  a 
surveyor's  "total  station"  with  electronic  data  logging  for  surface  locations 
and  elevations.  A  remotely  operated  vehicle  with  an  upward-looking  sonar,  a 
Mesotech  profiling  sonar,  and  a  thermal  drill  were  used  for  thickness  and 
sub-surface  details.  Table  1  gives  a  summary  of  the  data  collected  at  the 
seven  sites  documented  during  the  deployment.  A  detailed  description  of  the 
measured  data  as  well  as  a  description  of  the  measurement  techniques  can  be 
found  in  reference  18. 

On  Sunday,  September  29,1985,  the  participants  for  the  first  (one  week) 
leg  of  the  program  arrived  onboard  the  USCGC  POLAR  SEA  off  Barrow,  Alaska. 

The  POLAR  SEA  proceeded  approximately  40  nautical  miles  to  the  north  to  Site 
1  and  Site  2  (71.41  N,  157.03  W).  The  following  two  days  were  spent  preparing 
the  loads  data  acquisition  system  for  operation  and  conducting  on-ice 
measurements.  Ice  conditions  in  the  area  were  composed  of  large  first  year 
ice  floes  with  imbedded  multiyear  fragments.  The  heaviest  parts  of  the  floes 
ranged  from  5  to  8  meters  in  thickness.  A  total  of  eight  trial  rams  were  con¬ 
ducted  in  the  thinner  parts  of  these  floes  to  verify  the  correct  operation  of 
the  data  collection  system  during  the  evening  of  October  2.  This  was  followed 
by  two  ramming  tests  against  the  heaviest  ridge  sections.  These  rams 
completely  destroyed  both  of  the  principal  features  that  had  been  profiled. 

On  October  2,  a  helicopter  reconnaissance  flight  located  another  floe 
several  miles  away  that  contained  thicker  and  more  extensive  multiyear  ridges 
(Site  3  and  Site  4:  71.54  N,  157.22  W).  Ice  profiling  was  conducted  on 

October  3  and  October  4.  The  principal  multiyear  ridge  in  this  floe  was 
found  to  have  a  keel  of  approximately  42  ft  (13  m).  On  the  evening  of  October 
4,  this  ridge  and  several  adjacent  features  were  repeatedly  rammed  at  speeds 
up  to  nine  knots.  A  total  of  nine  rams  were  obtained  before  all  the  thick 
sections  of  ice  were  substantially  destroyed. 

A  personnel  transfer  was  conducted  off  Barrow  on  October  5.  Polar  Sea 
proceeded  approximately  60  nautical  miles  to  the  north  to  begin  the  second 
(eight  day)  leg  of  the  global  loads  program.  A  helicopter  reconnaissance  was 
conducted  while  the  ship  was  underway  and  Site  5  (72.23  N,  155.04  W)  was 
located.  Ice  in  the  vicinity  was  a  mixture  of  first  to  second  year  with  some 
small  multi-year  fragments.  Site  5  had  a  ridge  with  a  maximum  thickness  of 
approximately  46  ft  (14  m).  Two  ramming  tests  were  performed  on  October  6  at 
the  ice  edge,  however,  only  one  of  these  tests  was  recorded.  A  total  of  three 
rams  were  conducted  on  the  ridge,  on  October  8,  each  one  breaking  the  ridge  at 
various  locations.  Three  more  rams  were  conducted  at  surrounding  thick  ice 
features . 
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SITE 

LOCATION 

LATITUDE/ 

LONGITUOE 

SIZE  OF 
SURVEY 

AREA 

»  OF 

SURVEY 

POINTS 

MAX 

SAIL 

HEIGHT 

MAX 

KEEL 

DEPTH 

#  OF 

SONAR 

HOLES 

1 

71°41 ' NW 
157°03 '  W 

290 ‘ X290‘ 

172 

12.5 

36.1 

3 

2 

71°41’N 

157‘,03'W 

210'x270‘ 

140 

10.5 

28' 

0 

3 

7I°54‘ N 
157°22‘ W 

185'x215‘ 

160 

10.45 

30.1 

2 

4 

71<’54‘  N 
157°22‘W 

450‘x360‘ 

111 

14.4 

42.1 

2 

S 

72‘23‘N 

155°04'W 

440‘x270‘ 

185 

10.3 

35.7 

4 

6 

72“26'W 

154°56'W 

300'x200' 

79 

17.1 

37.5 

2 

7 

73°  N 
154°56‘U 

690‘x310‘ 

140 

12.7 

44.2 

3 

SITE 

#  OF 

THERMAL 

DRILL  HOLES 

TOTAL  DEPTH 
OF  THERMAL 

OR  ILL  HOLES 

ROV 

#  OF 
CORES 

TOTAL 

CORE 

LENGTH 

K/S 

1 

25 

470' 

0 

0 

2.67 

2 

71 

1074' 

2 

158“ 

2.67 

3 

62 

970' 

1 

174“ 

2.88 

4 

78 

1630' 

yes 

•  2 

336" 

2.92 

5 

100 

2161' 

yes 

3 

429“ 

3.47 

6 

11 

219* 

yes 

2 

90“ 

2.19 

7 

0 

0* 

yes 

2 

78“ 

3.23 

Table  1 

SUMMARY  OF  PROFILED  ENVIRONMENTAL  CONDITIONS 
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On  the  night  of  October  8,  a  large  ridge  at  Site  7  (73.00  N,  154.56  W) 
was  found.  In  the  following  two  days,  the  ridge  was  surveyed  and  found  to 
have  a  maximum  thickness  of  approximately  58  ft  (18  m).  The  floe  had  melt 
ponds  throughout  and  was  composed  of  first  and  second  year  ice.  A  total  of 
four  rams  were  conducted  on  the  ridge  on  October  10  each  one  breaking  the 
ridge  at  the  point  of  impact.  Fourteen  more  rams,  for  a  total  of  forty,  were 
conducted  on  ridged  ice  features  during  the  transit  back  to  Barrow,  Alaska, 
from  October  11  to  October  12. 

In  general,  the  typical  operation  was  to  locate  an  ice  feature,  posi¬ 
tion  the  ship  in  the  floe  next  to  the  feature  while  surveying  operations  were 
conducted,  and  then  to  ram  the  measured  feature  to  collect  the  ice  loading 
data  before  moving  on  to  the  next  site. 

It  was  anticipated  that  more  than  forty  ice  impacts  would  be  recorded 
during  the  deployment,  but  ice  conditions  were  not  as  severe  as  expected. 
Ideal  ice  conditions  for  the  tests  would  have  been  a  large,  thick  multiyear 
floe  that  could  be  rammed  repeatedly  without  breaking  apart.  Multiyear  ice 
was  only  encountered  in  small  floes  (some  with  a  reasonable  size  ridge  that 

was  surveyed)  that  were  imbedded  in  first  year  ice.  Much  of  the  old  ice  was 

only  second  year.  Table  2  summarizes  the  general  characteristics  of  each 
ramming  event. 

A  typical  ramming  test  consisted  of  the  following  sequence  of  events. 
The  icebreaker  cleared  a  path  to  ram  the  ice  feature.  The  vessel  was  moved 
perpendicular  to  the  ridge  feature  and  several  ship  lengths  away  before 

accelerating  towards  the  ridge.  Approximately  five  seconds  prior  to  impact 

the  data  acquisition  was  started  and  data  acquired  for  25  seconds.  Data 
collection  was  triggered  by  an  operator  viewing  a  video  display  of  the  bow  of 
the  vessel  and  the  immediate  area  ahead  of  the  ship.  The  measured  data  were 
transferred  to  a  floppy  disk  that  took  approximately  60  seconds  to  store. 

When  time  permitted  between  ramming  tests,  the  data  were  analyzed  and 
plotted. 
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TABLE  2 


OBSERVATIONS  OF  RAMMING  TESTS 


Magni tude. 

of 

Descri pti on 

Ram 

Measured 

of  Ice 

No. 

bate 

Symmetri c 

Load* 

Feature 

Comments** 

1 

10/02 

Y 

L 

Site  1 

2 

10/02 

Y 

L 

Site  2 

Good  Hit 

3 

10/04 

Y 

L 

Site  3 

12  m  Ridge 

4 

10/04 

Y 

L 

Site  3 

5 

10/04 

Port 

L 

Site  4 

Good  Heave 

6 

10/04 

Port 

S 

Site  4 

Roll  &  Pitch 

7 

10/04 

Y 

S 

Small  Floe 

Good  Acc.  ch 

8 

10/04 

Y 

S 

Small  Floe 

Floe  Submerged 

9 

10/04 

Y 

S 

10 

10/04 

Y 

S 

Small  Floe 

11 

10/04 

Y 

S 

Small  Floe 

12 

10/04 

Not  Recorded 

13 

10/06 

14 

10/07 

Y 

S 

Site  5 

Beaching  Force 

15 

10/07 

Y 

S 

Site  5 

16 

10/07 

Y 

L 

Site  5 

2nd  Hit  Hard 

17 

10/08 

Y 

L 

Site  6 

3  Hits 

18 

10/08 

Y 

L 

Site  6 

Hit  Ice  Knife 

19 

10/08 

Stbd. 

S 

Site  6 

Bounced  to  Port 

20 

10/08 

Y 

S 

Small  Floe 

21 

10/08 

Y 

L 

Small  Ridge 

22 

10/08 

Y 

L 

Small  Ridge 

23 

10/08 

Y 

L 

Cent.  Site  7 

3  Hinge  Break 

24 

10/10 

Y 

L 

Right  Site  7 

Floe  Submerged 

25 

10/10 

Y 

VL 

Left  Site  7 

Floe  Submerged 

26 

10/10 

Y 

VL 

Right  Site  7 

Floe  Submerged 

27 

10/11 

Y 

S 

In  Transit 

Small  Ridge 

28 

10/11 

Y 

S 

In  Transit 

Small  Ridge 

29 

10/11 

Y 

S 

In  Transit 

Small  Ridge 

30 

10/11 

Y 

S 

In  Transit 

Small  Ridge,  Ridge  Not  Broken 

31 

10/11 

Y 

M 

In  Transit 

Small  Ridge 

32 

10/11 

Stbd. 

M 

In  Transit 

Glanced  to  Port 

33 

10/11 

Y 

M 

In  Transit 

Large  Heave 

34 

10/11 

Y 

VL 

In  Transit 

Good  Hit 

35 

10/11 

Y 

M 

In  Transit 

36 

10/11 

Y 

S 

In  Transit 

Light  Hit 

37 

10/12 

Y 

S 

In  Transit 

Beaching  Force,  Not  Broken 

38 

10/12 

Y 

S 

In  Transit 

Impact  into  previous  bow 
print,  second  ram,  feature 

not  broken 

39 

10/12 

Y 

VL 

In  Transit 

Good  Impact,  third  ram. 

Ridge  broken 

40 

10/12 

Port 

L 

In  Transit 

Hit  Port  Side 

*  s. 

M,  L, 

&  VL  indicates 

small,  medium,  large,  and 

very  large  bow  force  loads. 

**  All 

rams 

resulted  in  failure  of  the 

ice  feature  except  where  noted. 
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4.  ANALYSIS  PROCEDURES  AND  RESULTS 


The  procedures  used  to  analyze  the  data  from  each  ramming  event  are 
summarized  here.  Appendix  B  gives  a  more  detailed  description  and  derivations 
of  the  various  equations  used. 

The  analysis  software  was  separate  from  the  data  acquisition  software. 
This  allowed  flexibility  during  the  data  collection  process  since  several  good 
rams  could  occur  a  few  minutes  apart.  The  separation  of  these  functions  (data 
acquisition  and  analysis)  meant  that  the  information  could  be  collected  and 
stored  for  future  analysis  without  missing  any  opportunities  to  collect  data 
during  sequential  rams. 

The  analysis  software  was  derived  from  the  program  written  for  the 
KIGORIAK  and  ROBERT  LEMEUR  impact  tests  conducted  in  1983.  Its  main  function 
was  to  calculate  and  plot  the  vertical  bow  force  time-history  acting  on  the 
POLAR  SEA,  and  determine  the  time  of  the  maximum  bow  force  together  with  the 
location  along  the  stem.  It  was  also  used  to  graph  the  shear  force  and  bend¬ 
ing  moment  distributions  at  any  time-step  during  the  25  second  sampling  inter¬ 
val.  In  addition  to  these  functions,  the  analysis  software  performed  a  number 
of  secondary  calculations  such  as  plotting  the  strain  time-history  from  any  of 
the  gages,  or  finding  the  location  of  the  neutral  axis  at  frames  instrumented 
on  two  separate  levels.  Appendix  C  contains  a  summary  of  the  program's  fea¬ 
tures  and  a  flow  chart  showing  the  branching  structure. 

During  the  analysis,  the  deck  bending  strain  time-history  at  each  gage 
location  is  calculated  and  plotted  for  every  ram.  Zeros  on  all  channels  are 
determined  by  averaging  the  data  obtained  just  prior  to  the  impact  and  sub¬ 
tracting  them  from  the  subsequent  measurements.  Stresses  are  calculated  by 
multiplying  the  results  by  a  calibration  factor  and  the  elastic  modulus.  The 
gages  that  were  placed  parallel  to  the  side  shell  of  the  vessel  are  multiplied 
by  another  transformation  factor  to  arrive  at  the  bending  stress  parallel  to 
the  centerline.  Appendix  A  contains  a  listing  of  the  gage  calibration  factors 
used,  while  Appendix  B  gives  the  derivation  of  the  centerline  stress  transfor¬ 
mation  for  each  instrumented  frame. 

A  sample  plot  of  the  bending  stress  at  frame  39  (usually  the  location 
of  the  highest  bending  stress)  is  shown  in  Figure  7.  During  the  ramming 
tests,  the  maximum  bending  stress  was  typically  determined  within  a  minute  of 
completion  of  the  ram.  These  values  were  always  well  below  the  yield  stress 
of  the  01  Deck  which  is  45,600  psi  (310  MPa).  A  histogram  of  the  maximum 
bending  stress  is  given  in  Figure  8  which  shows  that  the  highest  bending 
stress  recorded  is  about  6500  psi  (45  MPa). 
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TIME  (sac) 


Figure  9  illustrates  excitation  of  the  first  mode  natural  frequency  of 
the  hull  girder.  This  figure  shows  the  stress  time-history  measured  by  the 
starboard  bending  gage  located  at  frame  128  near  amidships  on  the  01  Deck. 

The  measured  frequency  of  3.0  Hz  is  within  8%  of  a  previously  measured  value 
of  2.75  Hz  [2]  and  very  close  to  the  2.9  Hz  computed  by  the  finite  element 
model  constructed  by  ABS. 


In  order  to  compute  a  bending  moment  from  the  bending  stresses,  it  was 
assumed  that  when  the  USCGC  POLAR  SEA  impacts  a  heavy  ice  feature  it  responds 
similar  to  a  beam  for  bending  within  the  centerline  plane.  The  bending  moment 
at  each  instrumented  frame  was  then  calculated  using  the  bending  stresses  and 
structural  properties  of  the  vessel. 


M 


where  o 
e 
E 
I 
Y 


Bending  stress  =  e  •  E 

Strain  (parallel  to  the  centerline) 

Elastic  modulus 

Transverse  sectional  moment  of  inertia 
Distance  between  gage  and  neutral  axis 
or  the  distance  between  gage  pairs  on  the  same  frame 


Referring  again  to  Figure  4  which  gives  the  locations  of  the  bending 
gages,  it  can  be  seen  that  several  of  the  frames  offer  a  couple  of  different 
methods  for  applying  this  formula.  At  cant  frame  44,  for  instance,  the 
stresses  at  the  gages  on  the  01  Deck  are  averaged  together  and  used  in  con¬ 
junction  with  either  the  gage  on  the  Third  Deck  or  the  First  Platform.  Or, 
for  that  matter,  the  01  Deck  gages  could  be  used  alone  along  with  their  ver¬ 
tical  distance  from  the  neutral  axis.  Generally  however,  a  gage  "couple"  was 
used  except  for  the  cases  where  bending  gages  were  installed  along  the  stem 
bar  and  a  few  of  the  other  centerline  gages.  These  particular  gages  were 
found  to  respond  to  the  local  load  of  the  ice  moving  down  the  stem  bar  or 
other  stress  concentration  influences  and  hence  were  not  used  in  the  calcula¬ 
tions. 


Once  the  bending  moment  distribution  along  the  length  of  the  ship  was 
obtained,  the  shear  force  was  computed  as  the  negative  slope  (derivative)  of 
the  longitudinal  bending  moment  curve.  Figure  3  shows,  generally,  how  these 
curves  appeared. 

Figure  3  also  shows  how  the  global  ice  force  is  related  to  the  shear 
diagram.  The  force  on  the  bow  was  calculated  by  the  addition  of  the  absolute 
value  of  the  greatest  shear  force  forward  and  aft  of  the  load.  The  location 
of  the  center  of  the  vertical  ice  load  was  estimated  from  the  measurements 
received  from  the  compression  gages  arranged  along  the  stem  bar.  At  any 
instant  in  time,  the  location  of  the  compression  gage  with  the  largest  com¬ 
pressive  strain  was  taken  as  the  ice  load's  location. 


16 


TIME  (sec) 


Figure  9 

BENDING  STRESS  TIME  HISTORY  NEAR  MIDSHIPS  (FRAME  128) 


This  entire  procedure  was  repeated  for  every  time  step  (0.01  seconds) 
for  the  duration  of  the  ramming  test  (25  seconds).  The  result  of  these  compu¬ 
tations  was  a  time-history  of  the  vertical  bow  force  during  the  ramming  event. 
Figure  10  gives  a  bow  force  time-history  plot  for  one  of  the  more  severe 
impacts.  Representative  rams  were  analyzed  onboard  the  vessel  using  prelimi¬ 
nary  estimates  for  the  sectional  inertias  and  locations  of  the  neutral  axes. 

It  was  anticipated  that  if  the  ice  load  occurred  forward  of  cant  frame 
17,  then  the  shear  force  would  be  estimated  by  the  multiplication  of  the  mea¬ 
sured  vertical  acceleration  and  the  mass  of  the  bow  section  forward  of  the 
load.  The  maximum  value  for  this  inertial  force,  however,  was  estimated  to  be 
approximately  30  LT  (0.3  MN)  which  is  less  than  the  uncertainty  expected  in 
computing  the  vertical  bow  force  and  was  therefore  neglected  (Section  6  dis¬ 
cusses  the  error  analysis). 

The  onboard  data  analysis  indicated  that  while  in  the  ramming  mode  the 
superstructure  of  the  POLAR  SEA  contributed  significantly  to  the  flexural 
stiffness  of  the  vessel.  This  was  apparent  when  the  calculated  bending  moment 
at  frame  55  (using  a  section  modulus  which  did  not  include  the  effect  of  the 
superstructure)  was  much  less  than  that  calculated  at  frame  39.  As  Figure  4 
shows,  frame  39  is  just  forward  of  the  superstructure  and  only  20.6  ft  (6.3  m) 
forward  of  frame  55.  The  bending  moment  distribution  for  this  portion  of  the 
ship  should  have  a  relatively  smooth  shape. 

The  calculated  bow  force  determined  from  the  discontinuity  in  the  shear 
curve  was  always  located  forward  of  the  superstructure  and  hence  unaffected  by 
the  sectional  properties  for  the  frames  under  the  superstructure.  Based  on 
this  observation,  it  was  decided  that  "effective"  sectional  properties  could 
be  found  for  these  frames  for  use  in  the  final  calculations.  The  location  of 
the  effective  neutral  axis  was  calculated  at  frames  where  the  bending  strain 
was  measured  at  two  levels  by  assuming  a  linear  stress  distribution  through 
the  cross-section.  The  point  where  this  distribution  passed  through  zero  was 
taken  to  be  the  effective  neutral  axis.  The  moment  of  inertia  for  each  of 
these  cross-sections  was  recalculated  from  the  ship's  drawings  based  upon  the 
new  location  for  the  neutral  axis.  With  the  assumption  that  the  ship  remains 
in  a  quasi-static  equilibrium,  the  areas  above  and  below  the  shear  diagram 
were  calculated  to  determine  if  they  were  equal.  This  was  done  as  a  check  on 
the  validity  of  the  recalculated  sectional  properties  since  they  were  not 
equal  using  the  original  estimates.  The  sectional  properties  for  the  frames 
including  the  superstructure  were  then  adjusted  to  bring  the  positive  and 
negative  areas  of  the  shear  graph  into  equilibrium. 
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Table  3  below  gives  a  listing  of  the  revised  neutral  axes  and  moment  of 
inertias  for  each  of  the  instrumented  frames  on  the  POLAR  SEA.  These  are  the 
values  used  in  the  final  calcul ations. 


TABLE  3 

SECTIONAL  PROPERTIES  FOR  THE  POLAR  SEA  (WAGB-11) 


LOCATION  NEUTRAL  AXIS  MOMENT  OF  INERTIA 

_ ft  (m) _ ft  \m  4) _ 


Frame 

128 

24.3 

(  7.41) 

11,586 

(100) 

Frame 

86 

25.8 

(  7.86) 

11,586 

(100) 

Frame 

55 

25.3 

(  7.7) 

15,062 

(130) 

Frame 

39 

28.2 

{  8.6) 

7,623.7 

(65.8) 

Cant 

Frame 

44 

34.8 

(10.6) 

3,058.8 

(26.4) 

Cant 

Frame 

35 

37.1 

(11.3) 

2,085.5 

(18.0) 

Cant 

Frame 

27 

40.4 

(12.3) 

1,181.8 

(10.2) 

Cant 

Frame 

17 

42.6 

(13.0) 

787.9 

(  6.8) 

As  a  check  on  the  validity  of  the  algorithm  used  to  calculate  the  total 
vertical  ice  force  on  the  bow,  a  comparison  was  made  with  the  beaching  force 
computed  from  hydrostatics.  A  couple  of  representative  rams  were  selected 
where  the  ship  came  to  a  complete  stop.  This  insured  that  the  longitudinal 
force  would  be  negligible.  At  the  instant  of  zero  headway,  the  vertical  bow 
force  was  computed  using  the  measurements  obtained  from  the  bending  gages. 

At  the  same  time  the  angle  of  trim,  <|>,  was  estimated  from  the  acceler¬ 
ometer  data.  The  hydrostatic  curves  were  then  used  to  obtain  the  moment  to 
trim  one  inch  (MTI)  and  the  beaching  force  calculated  from 

F  =  sin  $  •  LBP  •  MTI  /  d 

LBP  *  Length  between  perpendiculars 

where  d  is  the  approximate  distance  from  the  point  of  application  of  the  ice 
force  to  the  longitudinal  center  of  flotation.  Appendix  B  contains  further 
details  of  this  calculation.  Table  4  gives  some  of  these  results  from  which 
it  can  be  seen  that  the  values  obtained  using  the  different  methods  are  in 
fairly  good  agreement. 


20 


TABLE  4 


HYDROSTATIC  AND  CALCULATED  BOW  FORCE 


RAM 

NO. 


FORCE  FROM 
HYDROSTATIC  CURVES 
LT  (MN) 


FORCE  CALCULATED 
FROM  BENDING  GAGES 
LT  (MN) 


37  542  (5.4)  572  (5.7) 

39  371  (3.7)  361  (3.6) 


The  location  of  the  center  of  the  ice  force  (calculated  from  output  of 
the  stem  bar  compression  gages)  during  ramming  can  be  used  to  estimate  the 
impact  velocity.  A  sample  plot  of  the  calculated  location  of  the  load  versus 
time  is  shown  in  Figure  11.  The  slope  of  a  line  drawn  through  this  stepped 
curve  is  an  estimate  of  the  velocity  of  the  ice  movement  along  the  stem  bar. 
Correcting  for  the  angle  of  the  stem  bar,  an  approximate  val ue  for  the  ship 
impact  velocity  is  obtained.  A  comparison  between  the  impact  velocity  calcu¬ 
lated  from  the  location  of  the  load  time-history  and  the  velocity  measured 
from  the  doppler  speed  log  for  two  rams  is  shown  in  Table  5. 


TABLE  5 

COMPARISON  OF  IMPACT  VELOCITIES 

RAM  NO.  DOPPLER  SPEED  LOG  STEM  BAR  GAGES 

_ Loots _ (m/sec ) _ knots  (m/s  ec) 

8  7.0  (3.6)  8.6  (4.4) 

9  6.3  (3.5)  7.6  (3.9) 


The  difference  in  velocities  is  probably  due  to  the  nature  of  the 
ship-ice  interaction.  The  ice  moving  down  the  stem  bar  is  not  exactly  a  point 
load  and  does  experience  some  crushing  causing  the  point  of  maximum  loading  to 
shift  locations  within  the  ice  feature.  In  any  event,  the  velocities  given 
are  relatively  close. 

All  of  the  strain  data  was  analyzed  using  the  procedure  described  above. 
Table  6  summarizes  the  results  for  all  of  the  rams  and  gives  the  peak  vertical 
bow  force,  the  impact  velocity,  and  the  maximum  bending  stress  along  the  01 
Deck.  The  largest  bow  force  encountered  was  2506  LT  (24.97  MN)  during  ram 
number  39.  This  ram  also  obtained  the  highest  bending  stress  with  6078  psi 
(41.91  MPa)  in  compression  at  frame  39.  A  histogram  of  the  peak  vertical 
force  is  illustrated  in  Figure  12. 
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TABLE  6 


SUMMARY  OF  IMPACT  DATA  ANALYSIS 


RAM  IMPACT 


NUMBER 

VELOCITY 

kts  m/s 

MAX.  VERTICAL 

LT 

BOW  FORCE 

MN 

MAX.  BEND. 

psi 

STRESS 

MPa 

1 

8.3 

4.3 

851 

8.48 

-1530 

-10.55 

2 

8.9 

4.6 

1672 

16.66 

-3545 

-24.44 

3 

8.5 

4.4 

1516 

15.11 

-2515 

-17.34 

4 

5.9 

3.0 

552 

5.50 

-1008 

-  6.95 

5 

8.6 

4.4 

1402 

13.97 

-2235 

-15.41 

6 

5.9 

3.0 

804 

8.01 

-1748 

-12.05 

7 

6.1 

3.1 

1121 

11.17 

-1320 

-  9.10 

8 

7.0 

3.6 

806 

8.03 

-2054 

-14.16 

9 

7.0 

3.6 

1409 

14.04 

-3758 

-25.91 

10 

7.0 

3.6 

679 

6.77 

-1498 

-10 .33 

11 

8.6 

4.0 

717 

7.14 

-1697 

-11.70 

13 

4.0 

2.1 

528 

5.26 

-1043 

-  7.19 

14 

4.3 

2.2 

1610 

16.04 

-2554 

-17.61 

15 

5.7 

2.6 

1465 

14.60 

-1872 

-12.91 

17 

8.5 

4.4 

1624 

16.18 

-1997 

-13.77 

18(1) 

6.0 

3.1 

1473 

14.68 

-2312 

-15.94 

19 

4.3 

2.5 

1291 

12.86 

-2196 

-15.14 

20 

6.9 

3.5 

1151 

11.47 

-2055 

-14.17 

21 

8.5 

4.4 

1204 

12.36 

-1814 

-12.51 

22  f  1) 

6.0 

3.1 

1063 

10.59 

-1993 

-13.74 

23 

8.8 

4.5 

991 

9.87 

-1510 

-10.41 

24(1) 

7.0 

3.6 

841 

8.38 

-1523 

-10.50 

25 

7.5 

3.9 

1959 

19.52 

-4364 

-30.09 

26 

9.0 

4.6 

1463 

14.58 

-2716 

-18.73 

27(1) 

6.0 

3.1 

768 

7.65 

-1411 

-  9.73 

28(2) 

X 

919 

9.16 

-1516 

-10.45 

29(2) 

X 

613 

6.11 

-1357 

-  9.36 

3U(  2) 

X 

688 

6.86 

-1339 

-  9.58 

31 

5.0 

2.6 

1101 

10.97 

-1682 

-11.60 

32 

5.0 

2.6 

1009 

10.05 

-1524 

-10.51 

33 

11.5 

5.9 

2051 

20.44 

-2828 

-19.50 

34 

10.9 

5.6 

1418 

14.13 

-2445 

-16.86 

35 

6.6 

3.4 

679 

6.77 

-1352 

-  9.32 

36(1) 

9.0 

4.6 

934 

9.31 

-1846 

-12.73 

37 

6.0 

3.1 

650 

6.48 

-1433 

-  9.88 

38 

6.2 

3.2 

1923 

19.16 

-2483 

-17.12 

39 

8.5 

4.4 

2506 

24.97 

-6078 

-41.91 

40 

7.6 

3.9 

1603 

15.97 

-3192 

-22.01 

(1)  Data  from  bridge  velocity  indicator 

(2)  Data  not  available 
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5.  COMPARISON  OF  POLAR  SEA  RESULTS  WITH  PREVIOUS  REPORTS 


5.1  Peak  Vertical  Bow  Force  Versus  Impact  Velocity 

The  results  obtained  from  the  full-scale  impact  tests  onboard  the  USCGC 
POLAR  SEA  can  be  compared  with  previous  predictions  of  vertical  bow  force  on 
icebreaking  vessels.  Figure  13  gives  a  scatter  plot  of  the  vertical  bow  force 
versus  the  impact  velocity  for  all  impact  events.  Added  to  this  plot  are  sev¬ 
eral  maximum  vertical  bow  force  prediction  curves.  The  solid  curve  comes  from 
a  proposal  by  Johansson,  Keinonen,  and  Mercer  [14]  for  Arctic  Class  10  ves¬ 
sels.  They  felt  that  the  total  maximum  bow  force  was  largely  influenced  by 
the  ship's  speed  and  mass  and  gave  a  recommended  design  equation  of 

Fmax  =  V  -A0*9 

where  Fmax  =  maximum  force  in  MN 

V  =  ship's  speed  or  impact  velocity  in  m/s 
A  =  ship's  maximum  displacement  in 
millions  of  kilograms 

This  is  the  force  normal  to  the  hull,  and  so  the  vertical  bow  force 
component  would  oe  the  total  bow  force  times  the  cosine  of  the  angle  of  the 
stem  bar. 


^vert  =  ^max  *  cos  a 
0  9 

^vert  =  V  •  A  *  •  cos  a 

For  the  POLAR  SEA  the  displacement  is  close  to  11,000  LT  (11,170  MT)  at 
the  design  waterline  of  28  ft  (8.5  m)  and  the  stem  angle  5  feet  (1.5  m) 
below  this  waterline  is  about  22.5°.  The  corresponding  values  used  in  the 
above  equation  are 


A  =  11.17  millions  of  kg 
a  =  22.5° 

which  result  in  the  solid  curve  on  Figure  13.  This  curve  shows  Johansson's 
prediction  is  a  good  upper  bound  for  ramming  velocities  between  3.9  and  8.75 
knots  (2.0  and  4.5  m/s).  It  is  important  to  note  that  Johansson's  criteria 
was  intended  to  include  severe  ice  conditions  such  as  impacting  glacial  ice, 
while  the  ice  encountered  during  the  POLAR  SEA  trials  consisted  primarily 
of  multiyear  ridges  and  ridge  fragments  that  broke  upon  impact.  These  ice 
conditions  probably  account  for  the  lower  values  of  bow  force. 


A  second  comparison  can  be  made  with  the  full-scale  tests  conducted 
onboard  the  CANMAR  KIGQRIAK  in  1983.  The  inital  test  results  were  reported  by 
Ghoneim,  Johansson,  Smyth,  and  Grinstead  in  Reference  13.  They  developed  an 
envelope  curve  for  their  data  which  suggests  that  the  bow  force  is  propor¬ 
tional  to  the  square  root  of  the  impact  velocity. 

^vert  =  2 .34  /V  •  a  * 
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This  equation  is  indicated  by  the  dashed  line  in  Figure  13  with  the 
ship's  displacement,  A,  again  being  given  in  terms  of  millions  of  kg  to  give  a 
bow  force  in  terms  of  MN.  As  a  parametric  relationship  this  equation  was 
intended  by  Ghoneim,  et  al . ,  to  represent  only  a  possible  envelope  curve  based 
upon  the  portion  of  the  data  they  had  analyzed.  It  was  not  intended  to 
reflect  parameters  such  as  bow  shape.  The  ice  conditions  encountered  during 
the  KIGORIAK  tests  were  much  more  severe  than  the  POLAR  SEA  experienced  with 
many  of  the  rams  being  conducted  against  grounded  first  year  and  multiyear 
ridges.  As  the  graph  shows,  the  KIGORIAK  curve  certainly  does  provide  an 
envelope  to  the  POLAR  SEA  data  but  it  is  much  higher  than  Johansson's  curve. 
The  lower  vertical  bow  force  values  obtained  during  the  POLAR  SEA  tests  are 
again  probably  due  to  the  lighter  ice  conditions  although  it  would  be  diffi¬ 
cult  to  say  how  much  of  an  effect  the  different  bow  shapes  may  have  had. 

5.2  Vertical  Bow  Force  Time-Histories 


Typical  time-history  plots  of  the  vertical  bow  force  for  the  POLAR  SEA 
and  the  KIGORIAK  are  illustrated  in  Figures  14  and  15.  Ghoneim  and  Keinonen 
[6],  in  discussing  the  typical  ramming  scenario  for  the  KIGORIAK,  identify 
five  separate  phases.  These  are  the  approach  phase,  the  initial  impact  phase 
during  which  the  ice  crushes  and  the  bow  begins  to  ride  up  on  the  ice,  the 
slide  up  phase,  a  second  impact  phase  caused  by  the  knife  edge  contacting  the 
ice,  and  finally  the  slide  down  phase.  Figure  15  clearly  shows  the  two 
impact  phases  with  the  bow  force  dropping  from  2360  long  tons  (23.5  MN)  to 
zero,  then  rapidly  increasing  again  up  to  1200  long  tons  (12  MN) .  In  this 
case,  the  period  of  zero  ice  load  between  the  two  impacts  represents  the  bow 
rebounding  off  the  ice  surface  and  results  in  free  vibration  of  the  ship  until 
reloading  occurs. 

Figure  14  shows  a  typical  time-history  plot  from  the  POLAR  SEA  tests 
for  comparison.  After  the  initial  impact  of  around  1610  long  tons  (16  MN)  the 
bow  force  does  drop,  but  it  never  reaches  a  state  of  zero  ice  load.  That  is, 
the  bow-ice  contact  is  maintained  and  the  POLAR  SEA  does  not  "rebound"  as 
KIGORIAK  does.  In  fact  the  POLAR  SEA  bow  force  time-histories  do  not  show 
this  tendency  to  rebound  on  any  of  the  rams  analyzed  to  date.  The  displace¬ 
ment  of  the  POLAR  SEA  is  almost  1.7  times  that  of  the  KIGORIAK  which,  when 
coupled  with  a  different  bow  shape,  may  explain  the  difference  in  the  two 
types  of  response. 

Notice  that  in  the  plots  for  both  vessels,  a  relatively  constant  force, 
the  beaching  force,  is  achieved  at  the  end  of  the  ram.  This  is  the  force  that 
was  used  earlier  in  Table  4  for  the  comparison  of  hydrostatic  beaching  force. 
Also,  the  impact  duration  (length  of  time  of  the  first  peak)  is  approximately 
0.48  seconds  for  the  KIGORIAK  and  0.8  seconds  for  the  POLAR  SEA  for  these 
events.  Since  the  rams  occurred  at  different  velocities,  4.3  knots  (2.3  m/s) 
for  the  POLAR  SEA  and  9.5  knots  (4.9  m/s)  for  the  KIGORIAK,  a  more  extensive 
comparison  would  have  to  be  done  before  any  conclusion  could  be  reached  with 
regard  to  a  relationship  between  impact  duration  and  displacement. 
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Figure  15 

KIGORIAK  VERTICAL  BOW  FORCE  TIME  HISTORY 
(Ram  KR426  of  Reference  13) 


5.3  Longitudinal  Bending  Moment  and  Shear  Diagrams 

The  bending  moment  distributions  for  the  POLAR  SEA  and  the  KIGORIAK,  at 
the  time  of  maximum  force  for  a  large  impact,  are  shown  in  Figures  16  and  17 
respectively.  There  are  several  differences  between  these  graphs.  The 
location  of  the  load  on  the  POLAR  SEA  is  indicated  by  the  vertical  arrow  which 
intersects  the  bending  moment  curve  close  to  where  it  crosses  the  horizontal 
axis.  This  location  was  obtained  from  the  compression  gages  arrayed  along  the 
stem  bar  and  is  the  location  of  the  gage  reading  the  highest  compression.  For 
the  KIGORIAK  results,  a  "best  match"  procedure  was  employed  between  the 
bending  moment  and  shear  force  diagrams  to  estimate  the  center  of  the  load. 

The  maximum  bending  moment  of  the  POLAR  SEA  occurs  further  forward 
(approximately  75%  of  the  length  of  the  vessel  forward  of  the  stern)  than  the 
KIGORIAK's  (approximately  midships). 

The  shear  distributions  for  the  two  ships  (for  the  same  rams  used  in 
Figures  16  and  17  and  at  the  time  of  maximum  force)  are  shown  in  Figures  18 
and  19.  First  note  that  the  sign  convention  for  the  shear  force  is  opposite 
between  these  two  figures.  Since  bending  gages  were  installed  on  the  POLAR 
SEA  up  to  cant  frame  17,  which  was  forward  of  the  anticipated  maximum  load 
location,  the  bending  moment  and  shear  force  curves  could  be  calculated 
forward  of  the  load.  In  the  instrumentation  of  the  KIGORIAK,  however,  a 
slightly  different  approach  was  used  [13].  It  was  felt  that  since  the  bow 
force  was  concentrated  in  the  bow  area,  a  frame  instrumented  to  measure  the 
shear  force  just  aft  of  the  load  (frame  25  1/2  on  the  KIGORIAK)  would  be 
sufficient.  It  was  assumed  that  the  bending  moment  forward  of  the  load 
location  had  negligible  effect  on  the  computations.  Once  the  bending  moment 
and  shear  curves  were  obtained  up  to  frame  25  1/2  an  extrapolation  procedure 
was  used  to  obtain  the  bow  force  at  the  estimated  load  location.  It  is  for 
this  reason  that  Figure  19  does  not  show  any  shear  up  to  the  bow  (frame  30). 

Returning  to  the  shear  force  distribution  for  the  POLAR  SEA  (Figure 
18),  it  can  be  seen  that  around  the  location  of  the  load  the  shear  changes 
sign  over  approximately  50  ft  (15  m).  This  gives  a  rough  indication  of  the 
spreading  out  of  the  ice  load  over  the  extent  of  the  bow.  At  the  point  of 
maximum  vertical  force,  a  significant  amount  of  crushing  failure  has  occurred 
in  the  ice  feature  spreading  the  load  over  a  large  contact  area. 
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6.  EVALUATION  OF  THE  ACCURACY  OF  THE  GLOBAL  LOAD  MEASURING  SYSTEM 


There  are  a  number  of  potential  sources  of  error  that  should  be  con¬ 
sidered  in  order  to  estimate  the  overall  accuracy  in  computing  the  global  bow 
force.  Each  of  the  major  errors  is  investigated  in  turn  and  combined  with  the 
other  errors  using  the  standard  techniques  of  error  analysis  [19].  Starting 
with  the  equation  for  the  bending  moment  where  just  the  01  deck  bending  stress 
is  used  and  inserting  all  the  variables,  we  have 

M  =  a  *  I  =  e  •  E  •  I 
Y  Y 

M  =  e*  -  f  •  E  •  I 
Y 

where  f  is  the  transformation  factor  relating  the  strain  parallel  to  the  deck 
edge  to  the  principal  strain  along  the  centerline  and  Y  is  the  distance  from 
the  neutral  axis  to  the  gage  elevation  at  an  instrumented  frame.  The 
derivation  of  this  equation  is  contained  in  Appendix  B. 

The  error  associated  with  measuring  the  strain,  e' ,  can  be  estimated 
from  the  sampling  rate  and  the  accuracy  to  which  peak  amplitudes  of  a  signal 
are  measured.  Assuming  a  quarter  sine  wave  with  a  frequency  of  2.5  Hz  to 
represent  the  rise  in  strain,  as  mentioned  in  Section  2,  and  using  a  sampling 
rate  of  100  samples/second,  40  digital  samples  can  be  obtained  during  one 
cycle  at  the  highest  frequency.  The  digital  measurement  can  therefore  occur  a 
maximum  of  4.5°  away  from  the  peak  in  the  worst  case  (360°/(40x2) ).  This 
yields  a  maximum  error  in  sampling  the  peak  amplitude  of  ±0.31%. 

2 

Next,  the  expression  for  the  transformation  factor  f  contains  a  cos  0 
term  where  0  is  the  angle  the  strain  gage  is  positioned  off  of  the  center- 
line.  If  the  uncertainty  in  placing  the  gage  and  measuring  the  angle  is  about 
2°,  and  assuming  a  0  value  of  22°  (i.e.  near  the  bow),  then  the  uncertainty 
in  cos(e)  is  about  ±1.47%.  The  uncertainty  associated  with  the  transforma¬ 
tion  factor  would  then  be  twice  this  amount. 


Uncertainties  associated  with  the  moment  of  inertia,  I,  are  more  diffi¬ 
cult  to  determine.  For  the  frames  forward  of  the  superstructure  inertias  were 
computed  manually  from  the  ship's  plans  and  an  estimated  error  of  ±2.5%  was 
used.  The  neutral  axis  was  judged  to  be  accurate  to  within  ±0.5  ft.  Using  a 
value  of  40  feet  for  the  neutral  axis  near  the  bow,  then  the  resulting 
uncertainty  is  about  1.25%. 

Since  these  uncertainties  are  based  on  independent  measurements  they  can 
be  added  in  quadrature  to  arrive  at  the  uncertainty  in  calculating  the  bending 
moment . 

iM/M  =  [  (0.31)  2  +  (2x1.47)  2  +  (2.5)  2  +  (1.25)  2  ]0*5 

=  4.07% 
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Next,  ram  number  39,  the  maximum  one  recorded,  was  selected  to  observe 
how  the  uncertainty  in  the  bending  moment  propagated  through  the  equations  for 
shear  and  bow  force.  The  4.07  percent  computed  above  was  applied  to  the  four 
bending  moments  just  fore  and  aft  of  the  two  shear  forces  used  in  the 
calculation  of  the  vertical  bow  force.  The  uncertainty  for  just  one  of  these 
shear  forces  is  composed  of  the  errors  brought  about  by  the  uncertainty  in  the 
two  bending  moments.  These  two  errors  were  added  in  quadrature. 

4  S-j  -  [(A  Mi/4  X)  2  +  (4Mi+1/4  X)  2  ]°*S 

Finally,  the  uncertainty  for  the  two  shear  forces  were  also  added  in  quadra¬ 
ture  to  estimate  the  uncertainty  of  the  bow  force.  The  final  result  for  ram 
39  was  as  follows: 

4F  -  C(4  Si)2  +  (4Si+1)2]0-5 

4  F  -  246.8  LT 
F  =  2506  LT 

which  implies  4  F/F  «  9.8% 

Several  other  rams  were  analyzed  using  the  same  procedure  and  in  each 
case  the  uncertainty  was  less  than  10%.  This  overall  uncertainty,  however, 
takes  into  account  only  the  errors  associated  with  the  measurements  of  the 
individual  terms  that  make  up  the  expression  for  the  bending  moment.  Thus  the 
uncertainty  in  the  measurement  of  the  bending  moment  applies  only  to  values  of 
the  bending  moment  at  the  instrumented  frames.  It  does  not  include  any  error 
which  may  arise  from  measuring  the  bending  moment  at  a  finite  number  of 
points.  Therefore  an  additional  uncertainty  is  present  when  the  shear  force 
distribution  is  represented  by  the  slope  of  the  straight  line  segmented  bending 
moment  distribution.  A  more  reasonable,  but  qualitative,  assessment  of  the 
overall  uncertainty  in  the  bow  force  would  probably  be  ±15%. 

In  Section  4  the  vertical  bow  force  for  two  rams  was  compared  to  the 
force  obtained  from  a  hydrostatic  calculation.  The  larger  error  between  the 
two  cases  was  only  of  the  order  of  5%  which  lends  credibility  to  the  error 
estimates  computed  above. 

The  error  associated  with  the  location  of  the  center  of  the  ice  load  is 
unrelated  to  the  uncertainty  in  estimating  the  bow  force.  This  was 
determined  from  the  compression  gages  Installed  from  cant  frames  14  to  38. 

The  spacing  of  the  gages  allowed  the  load  center  to  be  estimated  with  an 
accuracy  of  ±  16  in  (±41  cm). 
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CONCLUSIONS 


The  following  is  a  list  of  the  conclusions  from  this  work: 

1.  The  use  of  bending  gages  provided  a  good  estimate  of  the 
longitudinal  bending  moment  and  shear  force  distributions.  This 
instrumentation  method  uses  fewer  strain  gages  than  attempting  to 
measure  the  shear  force  directly. 

2.  The  centerline  bulkhead  compression  gages  provided  an  excellent 
indication  of  the  location  of  the  center  of  the  ice  load.  They 
can  also  be  used  to  arrive  at  an  estimate  for  the  impact  velocity. 

3.  Global  impact  ice  forces  on  the  POLAR  SEA  are  not  localized  loads 
but  spread  over  much  of  the  bow. 

4.  The  superstructure  is  effective  in  bending  and  should  not  be 
ignored  in  design  calculations. 

5.  The  maximum  bending  stress  is  between  the  longitudinal  location  of 
the  forward  end  of  the  superstructure  and  the  bottom  of  the  stem. 

6.  The  loading  rate  was  measured  to  be  as  high  as  5000  LT/s  (50  MN/s), 
considerably  less  than  the  KIGORIAK's  loading  rate  of  15000  LT/s 
(150MN/S)  [13]. 

7.  The  vessel  did  not  "rebound"  after  the  first  impact  with  the  ice  as 
was  observed  in  the  KIGORIAK  tests  [13]. 

8.  The  duration  of  the  initial  impact  phase  of  a  ram  is  approximately 
0.8  seconds  for  a  4  knot  impact  velocity. 

9.  The  dominant  response  of  the  vessel  was  at  the  first  mode  of 
vibration  (3  Hz). 

10.  The  maximum  calculated  vertical  bow  force  was  2506  LT  (24.97  MN) . 

11.  The  maximum  bending  stress  measured  was  6078  psi  (41.91  MPa)  jn 
compression  at  frame  39  on  the  01  deck  which  is  well  below  the 
45,600  psi  yield  strength  for  the  deck  steel. 

12.  The  uncertainty  in  calculating  the  bending  moment  at  any  of  the 
instrumented  frames  is  approximately  ±4%.  The  propagation  of  this 
error  based  on  a  finite  number  of  points  results  in  a  bow  force 
uncertainty  of  ±  10%.  Since  the  bending  moment  distribution  should 
really  be  a  smooth  continuous  curve,  a  reasonable  estimate  for  the 
overall  uncertainty  is  more  likely  to  be  ±15%. 
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8.  RECOMMENDATIONS 


Recommendations  based  on  this  study  fall  into  several  categories; 
additional  analysis,  improvements  to  the  instrumentation,  and  additional  data 
collection. 

Additional  analysis  can  be  done  with  the  data  already  collected  to 
estimate  the  amount  of  local  ice  failure  during  the  ramming  event.  This 
report  sought  the  location  of  the  center  of  the  ice  load  from  the  stem  gage 
with  the  maximum  compressive  strain.  These  gages  could  be  analyzed  in 
conjunction  with  the  transverse  compression  gages  installed  on  cant  frame  35 
to  observe  how  the  contact  area  spreads  out  during  the  impact.  It  is  also 
recommended  that  further  analysis  be  done  on  the  unsymmetric  ramming  tests 
to  estimate  the  transverse  force. 

Future  instrumentation  programs  aboard  the  POLAR  Class  should,  wherever 
possible,  shift  centerline  strain  gages  outboard  to  the  deck  edge.  This  would 
reduce  the  problem  of  stress  concentration  influences  observed  with  gages 
installed  near  the  centerline  bulkhead.  The  data  from  the  two  separated  gages 
would  be  averaged  together  just  as  in  the  case  of  the  01  deck  gages.  It  would 
also  be  desirable  to  include  additional  instrumented  frames  in  the  bow  area  of 
the  ship  to  improve  the  definition  of  the  shear  fdrce  distribution  over  the 
region  of  the  ice  load. 

Additional  multiyear  ice  data  should  be  collected  with  the  POLAR  Class 
in  order  to  build  up  a  data  base  for  more  complete  analysis.  In  particular, 
ramming  events  collected  against  larger  ice  features  w  uld  help  to  validate 
the  possible  parameters  to  be  included  in  any  equations  describing  the  maximum 
bow  force  versus  impact  velocity. 
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APPENDIX  A 


SENSOR/CHANNEL  SPECIFICATIONS  AND  LOCATIONS 


This  appendix  gives  some  of  the  specific  details  associated  with  the 
operation  and  installation  of  the  data  acquisition  system  onboard  the  POLAR 
SEA. 


In  the  initial  configuration  of  the  system  maximum  expected  values  had 
to  be  estimated  for  the  strains  and  other  engineering  quantities  in  order  to 
calibrate  and  interpret  the  measurements  received.  Using  the  strain  gages  as 
an  example,  it  was  estimated  that  the  maximum  peak  strain  expected  to  occur 
with  the  highest  impact  loads  was  approximately  500  microstrain.  This  full 
scale  value  was  used  to  adjust  the  the  gain  setting  on  the  signal  filter/ 
amplifiers  to  provide  a  maximum  output  of  +10  volts  to  a  analog-to-digital 
subsystem.  The  A-to-D  system  used  was  a  HP  6940B  Multi  programmer  with  enough 
scanning  boards  to  handle  up  to  48  channels  of  data.  Associated  with  the 
analog-to-digital  process  was  a  conversion  factor  of  10  volts/2000  A-to-D 
units  which  places  some  limits  on  the  resolution  of  the  resulting  data.  For  a 
full  scale  of  497.512  microstrain/10  volts  on  all  the  strain  gages  this  reso¬ 
lution  corresponds  to  a  value  of  0.2487  microstrain/A-to-D  unit.  A  similar 
procedure  was  carried  out  for  the  accelerometers  and  the  doppler  velocity 
radar  to  arrive  at  the  calibration  factors  listed  in  Tables  A1  and  A2.  All 
of  the  strain  gages  were  wired  as  opposite-active  arm  bridges  and  had  an  exci¬ 
tation  of  10  volts. 

The  control,  storage  and  processing  data  collection  functions  were  per¬ 
formed  by  the  Hewlett-Packard  series  200  model  9816  desk  top  computer  with  an 
internal  memory  of  approximately  4.3  megabytes.  The  computer  was  used  to  pro¬ 
gram  and  receive  the  data  from  the  HP  6940B  Mul ti programmer.  Prior  to  each 
ramming  event  the  computer  programmed  the  HP  6940B  so  that  it  would  scan 
through  all  of  its  48  data  channels  (including  some  null  channels)  100  times 
each  second  for  25  seconds.  When  the  proper  cue  was  received  from  the  compu¬ 
ter  the  scanning  would  commence.  With  each  scan  the  multi  programmer  received 
the  filtered  data,  converted  it  to  digital  form,  and  transfered  it  to  the  com¬ 
puter's  memory.  Afterwards  the  measurements  were  transfered  from  computer 
memory  onto  a  3.5  inch  floppy  disk  for  storage;  one  ram  event  for  each  disk. 

The  signal  processing  filters,  analog-to-digital  multi  programmer,  and 
the  computer  were  all  set  up  in  the  windlass  room  aboard  the  POLAR  SEA.  This 
compartment  was  selected  because  of  its  central  location  with  respect  to 
laying  of  the  necessary  cables  to  each  of  the  strain  gages,  and  because  the 
out  of  the  way  location  would  not  interfere  with  shipboard  activities.  The 
equipment  rack  also  contained  a  video  screen  giving  the  view  from  aloft  conn 
of  the  bow  of  the  vessel  and  the  area  forward  of  the  bow.  Extra  components 
were  available  for  the  entire  system  in  case  the  failure  of  one  or  more  of  the 
essential  elements  occurred. 
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The  installation  of  electronic  equipment  and  strain  gages  onboard  the 
POLAR  SEA  was  carried  out  in  July  of  1985;  approximately  two  months  prior  to 
the  deployment.  After  arrival  in  the  Beaufort  Sea  and  before  any  ramming 
tests  were  conducted  the  entire  system  was  checked  out.  This  investigation 
led  to  four  strain  gages  that  were  found  to  be  defective  (probably  due  to 
water  exposure).  Three  of  the  gages  were  located  on  the  01  weather  deck  and 
replaced.  The  fourth  was  one  of  the  centerline  compression  gages  used  to 
locate  the  ice  load  along  the  stem  bar.  Unfortunately,  it  was  located  in  the 
lowest  part  of  the  3-E-O-W  tank  and  frozen  beneath  two  feet  of  ice.  It  was 
judged  that  ample  load  locating  capability  was  provided  by  the  rest  of  the 
compression  gages  and  so  no  attempt  was  made  to  replace  it.  The  refurbished 
gage  system  was  found  to  provide  clean,  noise  free  data  with  a  worst  case 
noise  of  less  than  20  millivolts  in  a  +10  volt  range. 
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TABLE  A1 


CHAN  £ 
1 


POLAR  SEA  IMPACT  TESTS  FALL  1985 


SENSOR  SPECIFICATIONS  FOR  RAMS  1  AND  2 


CHAN  ID 

CALIBRATION 

FACTOR 

GAGE 

FACTOR 

EXIT. 

VOLTAGE 

BRIDGE 

TYPE 

UNITS 

C-C24 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C26 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C28 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C30 

2.4 87E-6 

2.01 

10.00 

2.00 

strain 

C-C32 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C34 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C36 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C38 

2.487E-6 

2.01 

10.00 

2.00 

strain 

S-C27 

2.487E-6 

2.01 

10.00 

2.00 

strain 

S-C35 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B3-C35C 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B3-C44C 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B2-C17C  ' 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C20 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C18 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C16 

2.487E-6 

2.01 

10.00 

2.00 

strain 

C-C14 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B01-C17-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

BQI-C27-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B01-C35-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B01-C44-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B01-39-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B01-55-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B01-86-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B01-128-P 

2.487E-6 

2.01 

10.00 

2.00 

strain 

B3-39C 

2.487E-6 

2.01 

10.00 

2.00 

strain 

M-V  (speed) 

1.097E-2 

0.00 

6.09 

0.00 

knots 

M-A-T  (sway) 

6.225E-4 

0.00 

28.00 

0.00 

9 

M-A-L  (surge)  6.185E-4 

0.00 

28.00 

0.00 

9 

M-A-V  (heave) 

1.236E03 

0.00 

28.00 

0.00 

9 

C-C35-9P 

2.487E-7 

2.01 

10.00 

2.00 

strain 

C-C35-10S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

C-C35-1QP 

2.487E-7 

2.01 

10.00 

2.00 

strain 

C-C34-11S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

C-C35-11P 

2.487E-7 

2.01 

10.00 

2.00 

strain 

BIP-C44C 

2.487E-7 

2.01 

10.00 

2.00 

strain 

B01-C17-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

801-C27-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

801-C35-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

B01-C44-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

B01-39-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

BQ1-55-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

B01-86-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

B01-128-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

B3-55-P 

2.487E-7 

2.01 

10.00 

2.00 

strain 

B3-55-S 

2.487E-7 

2.01 

10.00 

2.00 

strain 

NULL 

NULL 
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TABLE  A2 

POLAR  SEA  IMPACT  TESTS  FALL  1985 


SENSOR  SPECIFICATIONS  FOR  RAMS  3  TO  40 


CHAN  # 

CHAN  ID 

CALIBRATION 

FACTOR 

GAGE 

FACTOR 

EXIT. 

VOLTAGE 

BRIDGE 

TYPE 

1 

C-C24 

2.487E-7 

2.01 

10.00 

2.00 

2 

C-C26 

2.487E-7 

2.01 

10.00 

2.00 

3 

C-C28 

2.487E-7 

2.01 

10.00 

2.00 

4 

C-C30 

2.487E-7 

2.01 

10.00 

2.00 

5 

C-C32 

2.487E-7 

2.01 

10.00 

2.00 

6 

C-C34 

2.487E-7 

2.01 

10.00 

2.00 

7 

C-C36 

2.487E-7 

2.01 

10.00 

2.00 

8 

C-C38 

2.487E-7 

2.01 

10.00 

2.00 

9 

S-C27 

2.487E-7 

2.01 

10.00 

2.00 

10 

S-C35 

2.487E-7 

2.01 

10.00 

2.00 

11 

B3-C35C 

2.487E-7 

2.01 

10.00 

2.00 

12 

B3-C44C 

2.487E-7 

2.01 

10.00 

2.00 

13 

B2-C17C 

2.487E-7 

2.01 

10.00 

2.00 

14 

C-C20 

2.487E-7 

2.01 

10.00 

2.00 

15 

C-C18 

2.487E-7 

2.01 

10.00 

2.00 

16 

C-C16 

2.487E-7 

2.01 

10.00 

2.00 

17 

C-C14 

2.487E-7 

2.01 

10.00 

2.00 

18 

B01-C17-P 

2.487E-7 

2.01 

10.00 

2.00 

19 

B01-C27-P 

2.487E-7 

2.01 

10.00 

2.00 

20 

B01-C35-P 

2.487E-7 

2.01 

10.00 

2.00 

21 

B01-C44-P 

2.487E-7 

2.01 

10.00 

2.00 

22 

B01-39-P 

2.487E-7 

2.01 

10.00 

2.00 

23 

B01-55-P 

2.487E-7 

2.01 

10.00 

2.00 

24 

B01-86-P 

2.487E-7 

2.01 

10.00 

2.00 

25 

B01-128-P 

2.487E-7 

2.01 

10.00 

2.00 

26 

B3-39C 

2.487E-7 

2.01 

10.00 

2.00 

27 

M-V 

1.097E-2 

0.00 

6.09 

0.00 

28 

M-A-T 

6.225E-4 

0.00 

28.00 

0.00 

29 

M-A-L 

6.185E-4 

0.00 

28.00 

0.00 

30 

M-A-V 

1.236E03 

0.00 

28.00 

0.00 

31 

C-C35-9P 

2.487E-7 

2.01 

10.00 

2.00 

32 

C-C35-10S 

2.487E-7 

2.01 

10.00 

2.00 

33 

C-C35-10P 

2.487E-7 

2.01 

10.00 

2.00 

34 

C-C34-11S 

2.487E-7 

2.01 

10.00 

2.00 

35 

C-C35-11P 

2.487E-7 

2.01 

10.00 

2.00 

36 

B1P-C44C 

2.487E-7 

2.01 

10.00 

2.00 

37 

BQ1-C17-S 

2.487E-7 

2.01 

10.00 

2.00 

38 

B01-C27-S 

2.487E-7 

2.01 

10.00 

2.00 

39 

B01-C35-S 

2.487E-7 

2.01 

10.00 

2.00 

40 

B01-C44-S 

2.487E-7 

2.01 

10.00 

2.00 

41 

B01-39-S 

2.487E-7 

2.01 

10.00 

2.00 

42 

B01-55-S 

2.487E-7 

2.01 

10.00 

2.00 

43 

B01-86-S 

2.487E-7 

2.01 

10.00 

2.00 

44 

B01-128  S 

2.487E-7 

2.01 

10.00 

2.00 

45 

B3-55-P 

2.487E-7 

2.01 

10.00 

2.00 

46 

B3-55-S 

2.487E-7 

2.01 

10.00 

2.00 

47  NULL 

48  NULL 
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UNITS 


strain 

strain 

strain 

strai n 

strai n 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

strain 

knots 

9 

9 

g 

strain 
strain 
strain 
strain 
strain 
strain 
strain 
strain 
strain 
strain 
strain 
strain 
strain 
strain 
strai n 
strain 


TABLE  A3 


POLAR  SEA  IMPACT  TESTS  FALL  1985 


SENSOR  LOCATIONS 

LONG  VERT 


CHAN# 

CHAN  ID 

ft 

(m) 

ft 

(m) 

1 

C-C24 

372.0 

113.4 

29.9 

9.1 

2 

C-C26 

362.5 

110.5 

29.2 

8.9 

3 

C-C28 

359.3 

109.5 

27.9 

8.5 

4 

C-C30 

356.3 

108.6 

26.9 

8.2 

5 

C-C32 

353.3 

107.7 

25.9 

7.9 

6 

C-C34 

350.3 

106.8 

24.6 

7.5 

7 

C-C36 

348.1 

106.1 

23.6 

7.2 

8 

C-C38 

344.5 

105.0 

22.0 

6.7 

9 

S-C27 

360.9 

110.0 

26.9 

8.2 

10 

S-C35 

348.8 

106.3 

22.6 

6.9 

11 

B3-C35C 

348.8 

106.3 

29.2 

8.9 

12 

B3-C44C 

335.6 

102.3 

29.2 

8.9 

13 

B2-C17C 

376.0 

114.6 

37.4 

11.4 

14 

C-C20 

371.4 

113.2 

30.8 

9.4 

15 

C-C13 

374.3 

114.1 

32.2 

9.8 

16 

C-C16 

377.3 

115.0 

32.2 

9.8 

17 

C-C14 

380.2 

115.9 

33.5 

10.2 

18 

801-C17-P 

376.0 

114.6 

59.4 

18 . 1 

19 

801-C27-P 

360.9 

110.0 

58.7 

17.9 

20 

501-C35-P 

343.3 

106.3 

58.4 

17. 8 

21 

BQ1-C44-P 

335.6 

102.3 

57.7 

17.6 

22 

B01-39-P 

316.9 

96.6 

56.8 

17.3 

23 

B01-55-P 

294.6 

89.8 

56.1 

17.1 

24 

B01-86-P 

253.6 

77.3 

54.8 

15.7 

25 

B01-128-P 

196.9 

60.0 

53.1 

16.2 

26 

83-39C 

316.9 

96.6 

28.5 

8.7 

27 

M-V 

193.6 

59.0 

52.5 

16.0 

28 

M-A-T 

369.1 

112.5 

49.5 

15.1 

29 

M-A-L 

369.1 

112.5 

49.5 

15.1 

30 

M-A-V 

369.1 

112.5 

49.5 

15.1 

31 

C-C35-9P 

348.8 

106.3 

25.9 

7.9 

32 

C-C35-10S 

348.8 

106.3 

24.6 

7.5 

33 

C-C35-10P 

348.3 

106.3 

24.6 

7.5 

34 

C-C34-11S 

348.8 

106.3 

23.3 

7.1 

35 

C-C35-11P 

348.8 

106.3 

23.3 

7.1 

36 

B1P-C44C 

335.6 

102.3 

21.0 

6.4 

37 

B01-C17-S 

376.0 

114.5 

59.4 

18.1 

38 

B01-C27-S 

360.9 

110.0 

58.7 

17.9 

39 

B01-C35-S 

348.8 

106.3 

58.4 

17.8 

40 

B01-C44-S 

335.6 

102.3 

57.7 

17.6 

41 

BQ1-39-S 

316.9 

96.6 

56.8 

17.3 

42 

BQ1-55-S 

294.6 

89.8 

56.1 

17.1 

43 

BQ1-86-S 

254.9 

77.7 

54.8 

16.7 

44 

801-128— S 

196.9 

60.0 

53.1 

16.2 

45 

83-55-P 

294.6 

89.8 

27.9 

8.5 

46 

B3-55-S 

294.6 

89.8 

27.9 

8.5 
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APPENDIX  B 


DETAILS  OF  THE  ANALYSIS  PROCEDURE 


CALCULATION  OF  BENDING  STRESS 

o(Rdg)  =  [(Data(Rdg,  Chn)  -  Base-line  (Chn)]x(Cal  ib(Chn )  )x (f  )x(E) 

o(Rdg)  =  8ending  Stress 

Rdg=  Reading  number  (100/second) 

Data  (Rdg,  Chn)  =  Oata  from  A/D  Subsystem 

Chn  =  Channel  number  (42  channels) 

Base-line  (Chn)  =  Base  line,  calculated  as  the  average  of 

100  readings  before  impact  occurs. 

Cal i b  (Chn)  =  Calibration  constant  (converts  A/D 

counts  to  strain) . 

E  =  Elastic  Modulus  30  x  106psi  (2.07E+11  Pa) 

f  *  Factor  for  transformation  of  stress 

parallel  to  side  shell  to  stress  parallel 
to  center  line  of  vessel. 

The  transformation  factor  f  which  relates  the  stress  parallel  to  the 
deck  edge  to  the  principal  stress  parallel  to  the  centerline  was  calculated 
from  the  equation  for  strain  at  a  point. 

2  2 

ex'  =  ex  •  (cos  9)  +  ey  •  (sine)  +  y  •  sine-cose 

If  ex  and  ey  are  the  strains  in  the  principal  directions,  then  the 
shearing  strain,  y,  is  equal  to  zero.  In  this  case  it  is  assumed  that  the 
principal  strain  ex  lies  along  or  is  very  close  to  the  centerline  of  the  ship. 

ex'  =  ex -(cos  g)  +  ey*(sine) 

ex'  =  ex«(coso)2  +  ey-(l  -  (cose)2) 

Using  ey  -  -  u«ex 

ex'  =  ex  •(cose)2  -  u*ex  •(  1  -  (cos  6)  2) 
ex'  =  ex  *C(COS0)  2  *(I  +  u)  -  u] 
ex  =  ex'/[(cos  9)  2  •(!  +  u)  -  u] 
ex  *  ex'  «f 
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where 


ex'  =  measured  strain  parallel  to  side  shell 

ex  =  principal  strain  parallel  to  the  centerline 

y  =  shear  strain  =  0 

0  =  angle  between  the  ship's  centerline  and  the 

tangent  line  to  side  shell  at  the  gage  location 
u  =  Poisson's  ratio  =  0.29 

Next  Hooke's  law  for  a  plane  stress  state  is  used  with  the  assumption 
that  for  the  ship  in  bending  the  stress  in  the  transverse  direction  is  much 
less  than  the  longitudinal  stress.  Therefore 

c(Rdg}  =  ex  -E 

o{  Rdg )  =  ex '  *f  *E 


CALCULATION  OF  BENDING  MOMENT 


Frame  55  M=( eZ- e3)x(E)x( I ) / D 

Frames  39,86,128  M=(  e01)x(E)x(  I) / Y 

Cant  frames  17,27,35,44 


M 

eOl 


bending  moment  (negative  if  01  deck  is  in 
compression) 

average  bending  strain  in  01  deck  at 
instrumented  frame 


e2 


average  bending  strain  in  2nd  deck  at 
instrumented  frame 


e3 


average  bending  strain  in  3rd  deck  at 
instrumented  frame 


E 


elastic  modulus  of  the  deck  plating 


I 


transverse  sectional  inertia 


Y 


distance  from  the  neutral  axis  to  the 
01  deck 


0 


vertical  distance  between  gages  located 
at  the  same  frame 


NOTE:  The  deck  bending  gages  located  on  the  centerline  were  not  used  in 

the  bending  moment  calculation,  because  local  stress  fields  occurred 
when  the  ice  force  was  located  at  the  frame  containing  these  gages. 
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CALCULATION  OF  SHEAR  FORCE  VALUES 


S(x) 

MI 

M2 

XI 


S(X)  =  -(M2-MI)/( X2-X1) 

the  shear  force  at  x=Xl  +  (X2-Xl)/2 

the  bending  moment  at  x=Xl 

the  bending  moment  at  x=X2 

the  longitudinal  location  of  bending 
moment  Ml 


X2  =  the  longitudinal  location  of  bending 

moment  M2 


LOCATION  OF  THE  CENTER  OF  THE  LOAD 

Loc  =  Loc_comp(Max_comp) 

The  location  of  the  center  of  the  ice  load  was  taken  to  be  at  the 
location  of  the  stem  bar  compression  gage  undergoing  the  largest  amount  of 
compression. 


CALCULATION  OF  VERTICAL  FORCE  ON  THE  BOW 
A.  If  bending  gages  forward  of  location  of  load 


Bow_f 

= 

ABS  CS(xl)] 

+  ABS  [S(x2)] 

Bow_f 

= 

the  vertical 

ice  force  on 

the  bow 

S(xl) 

= 

shear  force 

aft  of  the  location  of 

the 

S(x2) 

= 

shear  force 

forward  of  the 

location 

of 

ABS 

= 

absolute  value  of  quantity 

in  brackets 

load 

the  load 
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B.  If  no  bending  gages  forward  of  location  of  load 


Bow  f 

= 

ABS[S(xl) ]+(Mass )x(acc)  •  ABS(S(xl)) 

Bow  f 

= 

vertical  ice  force  on  bow 

S(xl) 

= 

shear  aft  of  location  of  load 

Mass 

= 

mass  of  bow  section  forward  of  ice  force 

acc 

- 

vertical  acceleration  of  bow  section 

NOTE:  The  maximum  inertial  force  [(mass )x (ACC ) ] ,  for  the  case  of  no  bending 
gages  forward  of  the  location  of  the  load,  is  approximately  30  LT 
(0.3  MN). 


CALCULATION  OF  BEACHING  FORCE  FROM  HYDROSTATICS 


Beach  f 


LBP 

MTI 


sin  4.  ♦  LBP  •  MTI/d 

beaching  force  on  the  bow  in  long  tons 
trim  angle 

length  between  perpendiculars  (inches) 

moment  to  trim  one  inch 

925  ft-LT  (at  a  draft  of  28  ft) 

distance  from  the  ice  force  to  the  longitudinal 
center  of  floation 


158  ft 


The  vessel's  angle  of  trim  when  all  headway  was  stopped  during  a  ramming 
event  was  estimated  from  the  accelerometer  measuring  the  longitudinal  accel¬ 
eration.  Before  the  ram  the  heave  accel erometer  indicates  a  Ig  acceleration 
while  the  sway  and  surge  accelerometers  indicate  zero  acceleration.  With  the 
vessel  in  the  beached  condition,  however,  the  heave  and  surge  accelerometers 
give  two  components  to  the  downward  lg  gravitational  acceleration. 

sin  <{>  ■  ax(measured  in  g's)/lg  =  ax 
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APPENDIX  C 


DESCRIPTION  OF  THE  GLOBAL  LOAD  ANALYSIS  SOFTWARE 


General  Description 

The  program  developed  for  the  USCGC  POLAR  SEA  Impact  Tests  was  derived  from 
the  program  written  for  the  M.V.  KIGORIAK  and  the  M.V.  ROBERT  LEMEUR  Impact 
Tests  conducted  in  1983.  The  main  function  of  the  program  is  to  calculate  and 
plot  the  global  vertical  load  time  history  on  the  bow  of  the  POLAR  SEA.  The 
following  calculations  can  be  performed  by  the  program: 

-  longitudinal  shear  force  and  bending  moment  distribution 

-  neutral  axis  location 

-  vertical  ice  force  on  bow 

-  cubic  spline  curve  fitting 

-  maximum  and  minimum  values  calculated 

-  longitudinal  location  of  maximum  and  minimum  bending  and  shear 

values 

-  location  of  center  of  ice  load  on  bow 

The  analysis  software  was  quite  versatile  in  its  mcde  of  application. 
While  a  full  length  report  containing  everything  from  graphs  of  the  output 
from  each  strain  gage  to  the  final  time-history  of  the  vertical  bow  force 

took  about  a  hour  to  generate,  the  program  allowed  for  much  quicker  results  if 

desired.  Once  the  data  disk  was  read  by  the  computer,  any  of  the  strain  gage 

time  histories  could  be  viewed  and  plotted.  Also,  any  time  segment  out  of  the 

25  second  sampling  interval  could  be  selected  for  calculation  of  the  bending 
moment,  shear,  or  bow  force  to  be  displayed  and/or  plotted.  This  latter  mode 
of  operation  reduced  the  computations  down  to  a  few  minutes  and  was  used  pri¬ 
marily  during  the  system  checkout  phase  to  determine  operational  readiness  of 
the  instrumentation. 

The  hardware  required  for  the  analysis  program  is  as  follows: 

-  HP  9816  Computer 

-  HP  9121  Disk  Drive 

-  HP  2225A  Think  Jet  Printer 

-  2.2  MBytes  of  Memory 

The  HP  Basic  3.0  operating  system  is  required  also. 

A  general  soft  key  menu  flow  chart  is  illustrated  in  Figure  Cl. 
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Figure  Cl 

GLOBAL  LOAD  ANALYSIS  FLOWCHART 


PROGRAM  USED:  CALC  I 

ANALYSIS  OATE/TIME:  20  Jan  1986/16:41:39 

ACQUISTION  OATE/TIME:  10  Oct  1985/18:01:55 


0-2 


CHNflSH  LTD 


1585  POLPR  SCPl  IMPACT  TT5TS  ( RRCTCC  CRNRBR  L^.  5 


TIME  S99 . 


PROGRAM  USED:  CALC! 

ANALYSIS  DATE/TIME:  Z0  Jan  1386/16:44:05 

ACQUISTION  OATE/TIME:  10  Oct  I  385/ 1 8 : 0 i : 55 
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PROGRAM  USED:  CALC! 

ANALYSIS  OATE/TIME :  Z0  Jan  1986/16:46:31 

ACQUISTION  OATE/TIME:  10  Oct  1985/18:01:55 
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1 305  polar  sen  iMPncr  rt5T3  iRRcrec  cnNftBn  LTC.  ; 


PROGRAM  USED:  CALC  I 

ANALYSIS  OATE/TIME:  20  Jan  1 S8S/ 1 6 : S I : 22 

ACOUISTION  OATE/TIME:  10  Oct  1985/18:01:55 
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t  sss  p<xhk  sen  irpnc:  reeTS  ( n scree  cflNnns  urc.  ? 


isas  polhr  sen  impact  tcsts  c nscrec  cnNncn  ltd.  j 


PROGRAM  USED:  CALC! 
ANALYSIS  OATE/TIME: 
ACQUISTION  OATE/TIME: 


20  Jan  I 98S/ 16:53:48 
10  Oct  1985/18:01 :S5 


0-12 


PROGRAM  USED:  CALC  1 

ANALYSIS  DATE/TIME;  20  Jan  I  986/ 16:55:01 
ACQUrsriQN  OATE/TIME:  10  Oct  1985/18:01:55 
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1385  P0U3R  SCR  IMPACT  TESTS  (RRCTEC  CHNHDR  LTC.  i 


PROGRAM  USED:  CALC1 

ANALYSIS  DATE/TIME:  20  Jan  1286/16:55:38 

ACOUISTION  OATE/TIME:  10  Oct  1985/18:01:55 
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1 985  POLFW  9Cfl  JtirmCT  TESTS  (IWCTEC  CANTO*  ltd.; 


1385  ROLPR  SEP!  TflAACT  TESTS  C RRCTEC  CPINRDR  LTD.  J 


PROGRAM  USED:  CALC1 

ANALYSIS  OATE/TIME:  20  Jan  1986/15:56:50 

AC0UI5TI0N  DATE/TIME:  10  Oct  1985/18:01:55 
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Isas  sen  inrncT  tests  (rt*CTCC  cnr^pn  ltd.  j 

RRM  +25 


l 

i 

\ 

i 


I 

i 

i 

i 


i 


isas  polrr  sen  impact  tests  t apctec  cnr-imm  lt3.j 


PROGRAM  USEO :  CALC! 
ANALYSIS  OATE/TIME: 
ACOUISTION  OATE/TIME: 


20  Jan  1586/16:58:03 
10  Oct  1285/18:01 :S5 
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1 305  POLAR  SEfl  IMPACT  TESTS  (  APCTEC  CANADA  LTD.  } 


l  305  POLAR  SEA  IMPACT  TESTS  l  AACTEC  CANADA  LTD.  1 


PROGRAM  USED:  CALC  1 

ANALYSIS  OATE/TIME:  20  Jan  1385/17:01:42 

ACOUISTION  DATE/TIME:  10  Oct  1985/18:01:55 


0-19 


1989  POLfW  5Cn  IMTflCT  TC5TS  (IWCTEC  CfVJfIDH  LTD.  3 


1989  POLFW  5CT  3T1PI9CT  TtS T5 


C  flffCTEC  CANADA  LTD.  3 


PROGRAM  USED:  CALC! 
ANALYSIS  OATE/TIME: 
ACOUISTION  OATE/TIME: 


20  Jan  1985/17:04:08 
10  Oct  1985/18:01 :5S 
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1333  POLAR  303  IMPACT  TESTS  (JWCTEC  CANADA  ltd.  3 
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TIME  S«5. 


1 333  POLAR  3C13  IMPACT  TESTS  (PMCTCC  CANADA  LTD.  3 


TIME  Svo . 


PROGRAM  USED :  CALC  I 
ANALYSIS  OATE/TIME: 
ACQUISTION  QATE/TIME 


20  Jan  1 986/ 17:07:45 
10  Oct  1985/18:01 -S5 


0-24 


